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ABSTRACT 

Classical lattice gas automata effectively simulate physical processes such as diffusion and 
fluid flow (in certain parameter regimes) despite their simplicity at the microscale. Mo- 
tivated by current interest in quantum computation we recently defined quantum lattice 
gas automata; in this paper we initiate a project to analyze which physical processes these 
models can effectively simulate. Studying the single particle sector of a one dimensional 
quantum lattice gas we find discrete analogues of plane waves and wave packets, and then 
investigate their behaviour in the presence of inhomogeneous potentials. 

PACS numbers: 03.65.-w, 02.70.-c, ll.55.Fv, 89.80.+h. 

KEY WORDS: quantum lattice gas; quantum cellular automaton; quantum computation. 



1 



Quantum mechanics of LGA I. 



David A. Meyer 



1. Introduction 

The first quantum lattice gas automaton (QLGA) appeared as Feynman's path integral for 
a relativistic particle in 1 + 1 dimensions [1] ; independently Riazanov constructed a 2 + 1 
dimensional QLGA as the path integral for the next higher dimensional Dirac equation [2] . 
In these formulations the quantum particle is conceptualized as evolving along spacetime 
trajectories, each of which is assigned a probability amplitude which is the product of a 
sequence of 'scattering' amplitudes describing the evolution of the particle during a single 
timestep. Thus these QLGA are discretizations of quantum mechanical processes. 

Feynman's path integral formulation of quantum mechanics reproduces the standard 
Schrodinger formulation of wave functions obeying partial differential equations [1] . These 
differential equations can be discretized directly, giving equations which Succi and Benzi 
naturally identify in the lattice gas paradigm as quantum lattice Boltzmann equations [3] . 
It is a familiar, although not often useful, observation that any numerical evolution of a 
discretized partial differential equation can be interpreted as the evolution of some cellular 
automaton (CA), if one allows the set of states to be R, or C, or for some very large 
N. Taking this perspective, Bialynicki-Birula constructs a model for quantum evolution — 
a linear unitary CA [4] — which is essentially equivalent to, although derived independently 
of, Succi and Benzi's equations. 

The equivalence of a QLGA simulation and the evolution of a set of quantum lattice 
Boltzmann equations/unitary CA depends on the equivalence of the path integral and 
standard formulations of quantum mechanics in the continuum. Our recent work explaining 
the necessity of non-unitarity in earlier attempts of Grossing, Zeilinger, et al. to construct 
homogeneous CA for quantum evolution [5] demonstrates this equivalence directly for 
the discrete models [6]. We also note that, in contrast to simulation with deterministic or 
probabilistic LGA, simulation with a QLGA requires evolution along all possible spacetime 
trajectories. This may be achieved (slowly) by evolution of the quantum lattice Boltzmann 
equation on a classical computer or, at present hypothetically (but rapidly), by simulation 
on a quantum computer. 

In fact, given the arguments that massive parallelism will optimize nanoscale quantum 
computer architecture [7] , it is plausible that the first useful quantum computation [8] will 
implement a QLGA simulation of some quantum mechanical process. This provides two 
reasons to pursue the project described in this series of papers: we want to explore not only 
quantum mechanical phenomena which can be simulated effectively by QLGA, but also how 
well, as Feynman suggested [9], a quantum computer might simulate physics. In addition, 
we expect the quantum mechanics of LGA to have implications for discrete models of 
fundamental physics: we have already found remarkable consequences of unitarity in linear 
[6,10] and nonlinear [11] QCA. 

We begin in the next section by recalling the model of [6] with which we will be 
working: the most general one dimensional homogeneous QLGA with a single particle of 
speed no more than 1 in lattice units. The local evolution rule for this model has two free 
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parameters: essentially the second measures the coupling between two copies of Feynman's 
original QLGA in which the first measures the 'mass' of the particle. 

This generalized QLGA is exactly solvable, just as is a single Feynman QLGA. In 
Section 3 we demonstrate this by finding the discrete analogues of plane waves in, and 
the dispersion relation for, our QLGA. We also show the results of simulations of the 
former — on a deterministic computer. 

We might imagine a one particle QLGA being simulated quantum mechanically by 
a ballistic electron in a solid state lattice [12] or as the 'low energy' sector of a line of 
dynamical quantum spins [9] ('low energy' meaning, e.g., the configurations with one spin 
up and the rest down). In the former case [13], and certainly if our interest is in the 
QLGA as a discrete approximation to the Dirac equation [6], it is natural to investigate 
wave packets representing a semi-classical quantum particle. We do so in Section 4. 

In Section 5 we show how to introduce an inhomogeneous potential into the model. 
Concentrating on finite square well potentials, we determine the dependence of the fre- 
quency/energy eigenvalues on the depth of the well and find that the eigenfunctions take 
the expected form. Finally, we utilize the results of Section 4 and show the results of 
simulations of a wave packet in a finite square well. 

We summarize our results in Section 6 and indicate the directions in which this re- 
search is continuing. 

2. The one particle QLGA 

A lattice gas automaton (LGA) should be envisioned as a collection of particles moving 
synchronously from vertex to vertex on a fixed graph (lattice) L: At the beginning of each 
timestep each particle is located at some vertex and is labelled with a 'velocity' indicating 
along which edge incident to that vertex it will move during the 'advection' half of the 
timestep. After moving along the designated edge to the next vertex, in the 'scattering' 
half of the timestep the particles at each vertex interact according to some rule which 
assigns new 'velocity' labels to each. For the purposes of this paper we will consider only 
one dimensional lattices L, isomorphic to the integer lattice Z or some periodic quotient 
thereof. In this case there are only two possible 'velocities': left and right. We will further 
restrict our attention to LGA with only a single particle; for some preliminary work on 
QLGA with multiple particles, see [6,14,15]. 

A QLGA is a LGA for which the time evolution is unitary. To make this precise 
we must first identify the Hilbert space of the theory. For a one particle QLGA in one 
dimension an ortho normal basis for the Hilbert space H is given by \x, a) (in the standard 
Dirac notation [16]), where x G L denotes position and a G {±1} denotes 'velocity'. At 
each time the state of the QLGA is described by a state vector in H: 




(2.1) 
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where the amplitudes ip a (t, x) E C and the norm of *&(t), as measured by the inner product 
on H, is: 

\ = ^$ a (t,x)%l) a (t,x). (2.2) 

x,a 

The state vector evolves unitarily, i.e., ^(t + 1) = U^S(t), where U is a unitary operator 
on H. Since the evolution is unitary, the inner product is preserved and (2.2) holds for all 
times if it holds for one; this allows the interpretation of ip a (t, x)ip a (t, x) as the probability 
that the particle be in the state \x, a) at time t [16,17]. As usual, therefore, the basis state 
vectors ^ = \x, a) correspond to 'classical' states — with probability 1 there is a single 
particle at x with 'velocity' a — and a generic state vector (2.1) is a superposition of these 
'classical' states, each of which has integer values (one 1, the rest 0) for the number of 
particles at each lattice site. (In general, the basis vectors for the n particle subspace of 
the QLGA Hilbert space are exactly the possible states of a classical deterministic LGA 
with n particles [15].) 

In order for the evolution to have the 'advection' interpretation described above, the 
basis vectors should evolve so that 

(x,a\U\y, 13)^0 (2.3) 

only when x = y + (3. This is equivalent to a condition on the amplitudes: 

ip a (t + l,x) = f a (ip-i(t,x + l),^ + i(t, x - 1)), 

where taking f a to be independent of x means that the QLGA is homogeneous in space. 
As the notation suggests, it is convenient to combine the left and right moving amplitudes 
at x into a two component complex vector ip(t, x) := (^_i(£, x), ip+i(t, x)) so that a state 
vector is written 

X 

We showed in [6] that the most general unitary evolution for a one dimensional QLGA 
with parity invariance* is unitarily equivalent to 

up to some overall phase which has no physical effect. Here the parameter 9 e K, or more 
precisely, tan 9, plays a role something like 'mass': when 9 = the particle travels only on 
the lightcone; as tan6> increases its probability for moving more slowly does also. 

Notice that just as in deterministic LGA in one dimension, the particle has a Z2 valued 
'Lagrangian' conserved quantity measuring the parity of its fiducial space coordinate. This 
'spurious' conserved quantity partitions the set of particles in a deterministic LGA into two 



I.e., invariance under x — > —x; also called reflection invariance. 
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decoupled gases [18] and in the QLGA denned by (2.4) it partitions the set of amplitudes 
ip(t, x) into two independent sets according to x + t (mod 2). This motivates consideration 
of the most general, no less local model which breaks this symmetry, namely 

ip(t + 1, x) = w-iip(t, x - 1) + w ip(t, x) + w +1 ip(t, x + 1), (2.5) 

where Wi G M2(C) are 2x2 complex matrices. In terms of the basis vectors \x,a), now 
(2.3) holds when x = y + (3orx = y, i.e., the particle can have nonzero amplitude to 
maintain its position. The condition that the global evolution, i.e., the matrix U, be 
unitary is expressed in terms of the Wi by the equations: 

w_ 1 w^_ 1 + w Wq + w +1 w+ 1 = I 

w w^_ 1 + w +1 Wq = (2-6) 
w +1 w i _ 1 = 0, 

together with their Hermitian conjugates [6]. Imposing also the condition of parity (re- 
flection) invariance on evolution of the form (2.5), we showed in [6] that the most general 
solution, up to unitary equivalence and an overall phase, is given by 

ism9\ ( cos# 0' 

w-i = cos p l n „ w+i = cosp . . a 

\ cos 9 \ % sin 6 , 

(2.7) 

sin 9 —i cos 9 
w = sin pi - a ■ Q 
\ —i cos 9 sin 9 

Here p G R is a coupling parameter breaking the 'spurious' symmetry. When p = 0, (2.5) 
reduces to (2.4), the QLGA which is unitarily equivalent to the models of Feynman [1], 
Succi and Benzi [3], and Bialynicki-Birula [4]. As tanp increases, the relative weight of wq 
increases and the particle has greater probability of maintaining its position, i.e., having 
zero velocity. This is the first indication of a symmetry between 9 and p which will become 
more explicit as we investigate the general QLGA of (2.5) and (2.7). 



3. Plane waves 



The local evolution rule (2.5) is linear so we expect the model to be exactly solvable. In 
[6] we solved the p = case by counting spacetime lattice paths in order to compute 
the propagator K a p(t, x; 0, 0) := (x, a|t/*|0, j3) explicitly. Lattice paths are more difficult 
to count when the particle has nonzero amplitudes for maintaining its position during 
each timestep. Avoiding this difficulty leads us to a more physical approach — finding the 
discrete analogue of plane waves in a QLGA. 

Recall that the QLGA is homogeneous, i.e., U commutes with the translation (shift) 
operator T defined by (Ti())(x) := ip(x + 1). Suppose L = Zjy. Then the eigenvalues 
of T are e lk for wave numbers k = 2-rm/N, n G {0, . . . , N — 1}, and the corresponding 
eigenvectors ty( k ' satisfy: 

ij {k) (x + l) = e lk ^ k) {x). (3.1) 
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Figure 1. Evolution of the n = 1 right moving Figure 2. Evolution of the n = 2 right moving 

plane wave on a periodic lattice with 9 = ir/3 plane wave on a periodic lattice with 6 = ir/3 

and p = 7r/4. and p = -it/ 4. 

Since [U,T] = and U is unitary, the are also eigenvectors for U with 

=e -*wk$(*) j (3.2) 

for some frequencies ujk £ K.. The eigenvectors are the discrete analogues of plane 
waves since they evolve simply by phase multiplication. 

Since the action of U is defined by (2.5), (3.1) and (3.2) imply that 

e- iLJh ^ k \x) = w-nlj ik) (x - 1) + w ^ k \x) + w +1 ^ k \x + 1) 

= (e'^w-! +w + e lk w +1 )vb {k \x) (3.3) 
=: D{k)^ k \x). 

Thus the e~ luJk are eigenvalues of D{k) e M 2 (C), i.e., solutions of 

det(D(k) -e~ iaJk I) = 0, (3.4) 

where / is the 2x2 identity matrix. Using the parametrization (2.7) of the u>i, (3.4) 
reduces to the condition 

cos u = cos k cos 9 cos p + sin 9 sin p. (3-5) 

For a given wave number k, (3.5) determines two frequencies ±u>k in terms of the rule 
parameters 9 and p. Call the corresponding eigenvectors of D(k) (normalized to have 
length 1/N) ^( fc ' ±1 ) (0) 6 C 2 , so that the corresponding plane waves are defined by (3.1) 
to be 

^(fc.e) . = J2^ (k,£) (0)e ikx \x). (3.6) 

X 

Figures 1 and 2 show the evolution of e = +1 (right moving) plane waves for n = 1,2. 
The probability ip*{t, x)ifj(t, x) (where ^{t,x) := t tjj(t,x)) of the particle being at x is 
constant in x (and t), so the vertical axis in the graphs shows the real part of ip-i(t,x). 
Even on such a small (N = 32) lattice this QLGA provides a very good approximation to 
continuum plane waves of long wavelength measured in lattice units. 
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Figure 3. The dispersion relation for 6 = n/3 Figure 4. The dispersion relation in the 'mass- 

and p = 7r/4. tt/12 < \u>\ < 5n/l2. less' case 6 = p = tt/6. \lu\ < 2n/3. 

Notice that when the wave number k increases, so does the frequency uj — the time 
period is shorter in Figure 2 than in Figure 1. Also the phase velocity ui/k decreases — the 
crest of the wave moves more slowly. In fact, (3.5) is the exact dispersion relation, giving 
the frequency in terms of the wave number. Figure 3 graphs the dispersion relation for 
the QLGA with the same rule parameters used in the simulations of Figures 1 and 2. The 
graph has reflection symmetry about both axes since (3.5) is invariant under both k — > —k 
and uj — > —uj. Each reflection alone changes the direction of the plane wave. When k = 0, 
uj = ±(6> — p); when k = ±7r, uj = ±(9 + p — ir). These values exemplify another symmetry 
of the dispersion relation — invariance under 9 < — > p; this is a symmetry in the QLGA rule 
space which is not realizable by a local unitary transformation. Figure 4 graphs the special 
case of equal rule parameters; here the dispersion relation passes through the origin. 

By comparison with plane waves in continuum quantum mechanics [16,17], we know 
that uj and k should be interpreted as being proportional to energy and momentum, 
respectively. Expanding the dispersion relation (3.5) around k = and uj = to second 
order, we find 

uj 2 = k 2 cos 9 cos p + 2 (1 - cos(# - p)) . (3.7) 
For a relativistic particle in the continuum, 

E 2 = p 2 c 2 + m 2 c 4 . (3.8) 

Comparing (3.7) and (3.8) suggests that the 1 — cos(6' — p) = case, i.e., the 9 = p case 
shown in Figure 4, corresponds to the particle being massless. 

Not only do the plane wave parameters ui and k bear the interpretation of proportion- 
ality to the conserved quantities energy and momentum, but they also label a complete 
set of (nonlocal) conserved quantities for the QLGA. Since T is orthogonal its eigenvectors 
i|r(fc,e) are orthogonal for distinct wave numbers k. Furthermore, D(k) is unitary, so its 
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eigenvectors tp( k,±1 >(0) are orthogonal for each k and hence so are the plane waves v[r( fe ' ±:L ). 
Since we normalized the eigenvectors of D(k) to have length 1/N, the plane waves (3.6) 
form an orthonormal basis for H which we denote by {|/c,e)}. Consider any state vector 
^ EH: 

x 

= ^yj(x) \k,e)(k,e\x) 
= ^2(^2(k,e\x)ij(x)^\k,e) 
The parenthesized expression is the amplitude of \k, e) in the new basis: 

X 

= E(E^ (fc ' e) (°) e ^^)V)^) 

x y 

= (^ (fc ' e) (0)) t ^V(^)e" lfcx (3.9) 

X 

=: (V> (fc ' e) (0))V(fc), 

where i^(k) is the discrete Fourier transform of ifj(x). The plane waves \k,e) evolve by 

phase multiplication so the probabilities tfj e (k)ifj e (k) are left invariant by the evolution. 
Since any initial state vector \l/(0) can be expressed in the plane wave basis this way, the 
existence of these conserved quantities is equivalent to exact solvability for this model of 
a one particle QLGA. 



4. Wave packets 

The plane waves (3.6) provide a starting point for constructing wave packets with localized 
position and particularized momentum. Consider the right moving plane wave with wave 
number /c in the position basis: 

v|>Oo,+i) _ y~^(fco, +1 )(Q} e ikoX \x). 

X 

In this discrete (and periodic) situation the binomial distribution is a convenient substitute 
for a Gaussian distribution, so to localize the particle we multiply the amplitudes by 
appropriate binomial coefficients: Let 

*-( 2 ;) _1 E^ +1, (0)^( x _ ;c ; +s/2 )w, (4.1) 
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Figure 5. Evolution of the fcg = 7r/4 wave packet (4.1) with width s = 32 for rule parameters 
9 = n/3, p=7r/4. The probability peak moves from x = 31 at t = to x = 54 at t = 49; thus 
the group velocity is approximately 23/49 ~ 0.47. 

for even s < N, where the inverse binomial coefficient outside the sum is the requisite 
normalization factor. This wave packet is localized around Xq, having support on the 
interval [xq — s/2, xq + s/2]. Figure 5 shows the evolution for wave number ko — 7r/4 and 
width s = 32 on the lattice The rule parameters are the same as those used in the 
simulations shown in Figures 1 and 2. In contrast to those graphs, the vertical axis in 
Figure 5 shows the probability that the particle is in the state \x). 

This simulation shows that the ko = 7r/4 wave packet moves with well defined group 
velocity to the right. The result is just what we would expect by analogy with the contin- 
uum situation and can be analyzed in the same way, by transforming to the \k,e) basis. 
Using (3.9) we compute the amplitudes in this basis: 



$ 6 (fc) = (0)) T ( ^ ) ^ (fco,+1) (0)e lkox 



](2s 

s) ^ T \x — xq + s/2 



x 



e -ik 3 



(^ (fc ' e) (o))V (fco,+1) (o) 



U( fe o,+i)rn^ 2sN \ ( s \j(k -k)x 



S J ' \X — Xn + s/2 

' x x ' 

1 



e 



(V> (fc ' e) (0)) V (/c °' +1) (0) f J e i{ - ko ~ k){x °~ s/2 \l + e^o-Vy 
(V; (fc ' e) (0))V (/£o ' +1) (0)f 2S ') 1 e ^- fc )^2 s cos s (^— -) (4.2) 
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M 2 




Figure 6. Evolution of the fcp = 7r /4 wave packet (4.1) with width s = 8 for rule parameters 
9 = 7r/3, p = 7r/4. This wave packet disperses more rapidly than the one shown in Figure 5: 
The peak probability at the end of the simulation is less than half of the initial peak probability; 
left moving ripples carrying off some of the probability are also visible. 



The amplitudes (4.2) give probabilities peaked around k = ko, so this is also a wave packet 
in momentum space. As usual, the group velocity is the slope of the dispersion relation 
(3.5), i.e., du/dk\ ko , which is y/9-2y/E/4 fa 0.49 for the values used in the simulation of 
Figure 5; this is in good agreement with the measured value of approximately 0.47. 

The width of the peak in (4.2) depends inversely on s: as s decreases, i.e., the width 
of the wave packet in position space decreases, the width of the momentum peak increases. 
The simulation in Figure 6 shows the evolution of a wave packet with width s = 8. We 
note that while the group velocity is the same as in Figure 5, there is substantially more 
dispersion, indicating a greater interval of contributing wave numbers. This is a general 
result, not depending on the specific form of our wave packet; the reciprocity relation for 
the discrete Fourier transform has consequences similar to those of the uncertainty relation 
for the continuous Fourier transform [19]. 

Figure 7 shows a simulation of a wave packet built from the plane wave with smallest 
nonzero wave number on the lattice: ko = n/32. The horizontal tangent to the graph 
of the dispersion relation at k = 0, as shown in Figure 3, indicates that the group velocity 
of this wave packet will be small. Furthermore, even with width s = 32 the wave number 
interval includes the left going modes whose presence is visible in Figure 7; the consequence 
is an interference pattern and no very well defined group velocity. 



10 



Quantum mechanics of LGA I. 



David A. Meyer 




Figure 7. Evolution of the fcg = 7r/32 wave packet (4.1) with width s = 32; the rule parameters 
are still 8 = 7r/3, p = 7r/4. This wave packet disperses even more rapidly than the one shown 
in Figure 6: left moving waves carry off some of the probability and an interference pattern is 
created. 




Figure 8. Evolution of the fcg = 7r/32 wave packet (4.1) with width s = 32 for rule parameters 
6 = 7r/6 = p. This wave packet disperses very little and has group velocity close to 1 in lattice 
units. 
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Finally, Figure 8 shows the evolution of the same wave packet but for the rule pa- 
rameters 6 = 7r/6 = p whose dispersion relation is graphed in Figure 4. Here the group 
velocity is close to 1 in lattice units, even for ko as small as 7r/32; the particle is indeed 
'massless'. There is almost no dispersion in this simulation; the probability contained in 
left going modes is nonzero, but too small by several orders of magnitude to be visible in 
Figure 8. 

5. Potentials 

The one particle QLGA described in Section 2 is the most general homogeneous model for 
particle speeds no more than 1. To simulate physical systems (or to do useful computation), 
some inhomogeneity must be introduced. In each of the equations (2.6), which express the 
unitarity condition, all the Wi correspond to the scattering/interaction at a single lattice 
point, as do all the w\. In the first equation these are the same lattice point, while in the 
second and third they are different. Thus if Wi(x) = Wi, constants independent of x, solve 
these equations, so do e~ l ^ x ^Wi. As observed already by Feynman [1] and Riazanov [2], 
such an x dependent phase realizes an inhomogeneous potential in the continuum limit 
of the discrete path sum for the Dirac equation. Here we investigate its effects on the 
quantum mechanics of our LGA, expecting them to be similar to those in the continuum 
limit. 

For simplicity, we restrict our attention to a finite square well potential, i.e., 

( x / e-^Wi if N/4 < x < 3A/4; 
L Wi otherwise, 



where the w\ are defined by (2.7). We begin 
for </>. 

Recall that the frequency (or energy) 
eigenvalues u are doubly degenerate ex- 
cept for those with the largest and small- 
est absolute value. (See Figure 3, where 
each horizontal line intersecting the graph 
of the dispersion relation does so at two 
points except when tangent to the maxi- 
mum or minimum of either branch of the 
curve.) As with any perturbation to the 
evolution, we expect the introduction of 
an inhomogeneity in the potential to re- 
solve the degenerate eigenvalues. Figure 
9 shows that this is indeed the case: as <p 
increases away from the eigenvalues u of 
U increase and the degenerate ones split. 



by considering the effect of different values 



CO 







71 








7C/2 




-31 








V / // - K 





<p 



Figure 9. The eigenvalues u> of U for a square 
well of depth <fi and width N/2 on a lattice of size 
N = 8 with 6 = tt/3 and p = tt/4. 

8; Figure 10 shows 



The eigenvalues in Figure 9 have been computed for only = 
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Figure 10. The eigenvalues uj of U for a square Figure 11. The eigenvalues u> of U for a square 

well of depth <fi and width N/2 on a lattice of size well of depth <j> and width N/2 on a lattice of size 

N = 32 with 6 = n/3 and p = ir/4. N = 32 in the 'massless' case 8 = n/6 = p. 

the results for N = 32. On the larger lattice it is clear what happens: the horizontal 
bands of frequency /energy eigenvalues correspond to the eigenvalues of the unperturbed, 
homogeneous system, while the diagonal bands of eigenvalues correspond to the same ones, 
but shifted by the depth <p of the square well. The periodicity along the frequency axis 
shown in these graphs is a symptom of the ambiguity in the definition of energy due to 
discrete time evolution [20]. The graphs in Figures 9 and 10 have been computed for the 
QLGA with parameter values 9 = n/3, p = n/A : the dispersion relation for which is shown 
in Figure 3. Repeating the calculations for the 'massless' case, with dispersion relation 
shown in Figure 4, gives the frequency /energy eigenvalue plot shown in Figure 11. The 
degenerate levels still split, but much less than before for the same (f> values, and the part of 
the band structure resulting from the nonzero minimum positive frequency in the massive 
dispersion relation vanishes. 

Now consider the eigenvectors of U, namely the eigenfunctions for our discrete version 
of a finite square well. Since U is no longer translation invariant we do not have an 
equation like (3.3) to solve for the eigenfunctions analytically. Rather than developing 
a cross boundary matching method as is used in the continuum problem for a periodic 
square well potential [21], here we simply find the eigenvectors of U numerically. Figure 
12 shows the eigenfunctions corresponding to the three smallest positive eigenvalues for 
the QLGA with 9 = tt/3 and p = it /A. The depth of the square well is <6 = 7r/24. We see 
exactly the lowest modes we would expect from our experience with such a potential in 
the continuum. As the energy of the eigenfunction increases there is greater probability 
that the particle is outside the well — in the region of higher potential. Figure 13 shows 
an eigenfunction which is approximately a plane wave in both regions: it has larger wave 
number in the well than outside it. For analytic results on the closely related problem of 
a step potential, and some discussion of their consequences for the physical interpretation 
of QLGA, see [15]. 
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Figure 12. The three eigenfunctions of U with 
smallest positive eigenvalues: 7r/12 < 0.2622 < 
0.2634 < 0.2653, for a square well of depth tt/24 
and width N/2 on a lattice of size N = 256 with 
9 = 7t/3 and p = it/4. 
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Figure 13. An eigenfunction for a particle with 
eigenvalue 0.3985 < 57r/12 large enough not to 
be confined completely to the square well of the 
previous figure. Outside the square well the wave 
number decreases and the probability increases. 



Finally, suppose we prepare one of the semiclassical wave packets studied in Section 
4 in a finite square well. Using the dispersion relation (3.5), we find that the ko = it/ 4, 
width s = 32 wave packet (4.1) of Figure 5 has peak frequency u = Cos _1 ((l + \/6)/4) 
for rule parameters 9 = tt/3, p = ir/4. ujq is just a little larger than it/6 so we would not 
expect a square well of depth = tv/6 to contain this wave packet. Figure 14 shows a 
simulation of this situation on a lattice of size N = 64: the wave packet continues past the 
right edge of the square well at 3iV/4 with only a small amount of internal reflection. 

Increasing the depth of the square well should have the effect of increasing the amount 
of internal reflection of the wave packet. Simulations demonstrate that this is indeed the 
case. When the depth of the square well is <fi = tt/4, Figure 15 shows that the wave 
packet splits as it scatters off the right wall of the square well. With greater probability 
the particle is reflected back into the well, but it also has a substantial probability of 
continuing to the right. The wave packet which continues to the right does so at a reduced 
group velocity as we can see by the fact that the reflected wave packet travels back across 
the well, reaching the left wall at x = N/4 = 16 at the end of the simulation shown, before 
the transmitted wave packet travels the same distance rightwards. 

Finally, when the depth of the square well is increased to = 7r/3, Figure 16 shows 
that almost the entire wave packet is reflected back into the well by the right wall. In this 
case there is only a very small probability that the particle has sufficiently large energy to 
escape the well. 



6. Discussion 

Unitarity is a very restrictive constraint on the local scattering rule for a QLGA with a 
single particle of bounded speed. When the bound is 1 in lattice units, there is a two pa- 
rameter family of reflection invariant one dimensional local rules, given by (2.5) and (2.7). 
It is already remarkable that the Dirac equation arises as a continuum limit of this QLGA 
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Figure 14. Evolution of the fc = 7r/4 wave packet (4.1) with width s = 32 for rule parameters 
9 = 7r/3, p = 7r/4 in a square well of depth <p = n/6. There is very little reflection as the wave 
packet passes the right wall of the square well at x = 3N/4 = 48. 

when p = 0. In this paper we have begun to investigate the quantum mechanics of the gen- 
eral two parameter rule. We find that even without going to a continuum limit the QLGA 
reproduces the quantum mechanical phenomena of plane waves and wave packets obeying 
a dispersion relation (3.5). Furthermore, the model straightforwardly accommodates the 
inclusion of inhomogeneous potentials. The eigenvectors of the evolution matrix give the 
quantum mechanical eigenfunctions for the lattice gas particle, and simulations exhibit the 
semi-classical evolution of a wave packet, in the presence of a square well potential. 

Taking a QLGA seriously as a possible model for quantum computation by, for ex- 
ample, ballistic electrons in a lattice of solid state nanostructures, raises many additional 
questions, some of which will be addressed in subsequent papers in this series: Inhomo- 
geneity of the substrate can be incorporated in the model by varying the rule parameters 
while maintaining global unitarity. Finite, non-periodic, boundary conditions can be im- 
posed similarly [22]. Higher dimensional [2,4,10,14] and multiparticle [6,14,17] models can 
also be constructed. Decoherence is the crucial problem for quantum computers [8], par- 
ticularly in the solid state [23] . QLGA provide an extremely convenient arena in which to 
model this problem [24]. Finally, the question of for which quantum computational tasks 
QLGA are best suited deserves serious investigation. 
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(figure available from author) 



Figure 15. Evolution of the same wave packet with the same rule parameters as in Figure 14, 
but now in a square well of depth <f> = it /A. There is both reflection and transmission as the 
wave packet scatters off the right wall of the square well. 



(figure available from author) 



Figure 16. Evolution of the same wave packet with the same rule parameters as in Figures 14 
and 15, but now in a square well of depth 4> = n/3. This well is deep enough that the wave 
packet is almost entirely reflected by the right wall of the square well. 
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